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Abstract We present a Lagrange-Galerkin scheme free from numerical quadrature 
for the Navier-Stokes equations. Our idea is to use a locally linearized velocity and 
the backward Euler method in finding the position of fluid particle at the previous 
time step. Since the scheme can be implemented exactly as it is, the theoretical stabil¬ 
ity and convergence results are assured. While the conventional Lagrange-Galerkin 
schemes may encounter the instability caused by numerical quadrature errors, the 
present scheme is genuinely stable. For the P 2 /P 1 - and Pi+/Pi-finite elements opti¬ 
mal error estimates are proved in (II 1 ) x l 2 (l 2 ) norm for the velocity and pressure. 
We present some numerical results, which reflect these estimates and also show the 
genuine stability of the scheme. 

Keywords Lagrange-Galerkin scheme • Finite element method • Navier-Stokes 
equations • Exact integration 

Mathematics Subject Classification (2000) 65M12 • 65M25 • 65M60 • 76D05 • 
76M10 


1 Introduction 

The purpose of this paper is to present a Lagrange-Galerkin scheme free from nu¬ 
merical quadrature for the Navier-Stokes equations and to prove the convergence. 
The Lagrange-Galerkin method, which is also called characteristics finite element 
method or Galerkin-characteristics method, is a powerful numerical method for flow 
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problems, having such advantages that it is robust for convection-dominated prob¬ 
lems and that the resultant matrix to be solved is symmetric. It has, however, a draw¬ 
back that it may lose the stability when numerical quadrature is employed to integrate 
composite function terms that characterize the method. Our scheme presented here 
overcomes this drawback. 

Lagrange-Galerkin schemes for the Navier-Stokes equations have been devel¬ 
oped in mmiMTiuiriW7iri9i : see also bibliography therein. After convergence 
analysis was done successfully by Pironneau ED in a suboptimal rate, the opti¬ 
mal convergence result was obtained by Stili |T9l . Optimal convergence results by 
Lagrange-Galerkin schemes were extended to the multi-step method by Boukir et al. 
0 , to the projection method by Achdou-Guermond |Tj and to the pressure-stabilized 
method by Notsu-Tabata in. All these results of the stability and convergence are 
proved under the condition that the integration of the composite function terms is 
computed exactly. Since it is difficult to perform the exact integration in real prob¬ 
lems, numerical quadrature is usually employed. It is, however, reported that insta¬ 
bility may occur caused by numerical quadrature error for convection-diffusion prob¬ 
lems in rnnmm . We observe such instability occurs for the Navier-Stokes 
equations by numerical examples in this paper. 

Several methods have been studied to avoid the instability in l4l ll2l[T6llT7l[23l . 
The map of a fluid particle from the present position to the position a time increment 
At before (the position is often called foot along the trajectory) is simplified. To find 
the foot of a particle is nothing but to solve a system of ordinary differential equations 
(ODEs). Morton et al. |T2] solved the ODEs only at the centroids of the elements, 
and Priestley iei did only at the vertices of the elements. The map of the other 
points is approximated by linear interpolation of those values. It becomes possible 
to perform the exact integration of the composite function terms with the simplified 
map. Bermejo et al. 0 used the same simplified map as im to employ a numerical 
quadrature of high accuracy to the composite function terms for the Navier-Stokes 
equations. Tanaka et al. 11231 and Tabata-Uchiumi l22l approximated the map by a 
locally linearized velocity and the backward Euler approximation to solve the ODEs 
for convection-diffusion problems. The approximate map makes possible the exact 
integration of the composite function terms. 

In this paper we prove the convergence of a Lagrange-Galerkin scheme with the 
same approximate map as H22II23I in the P 2 /P 1 - or Pi +/Pi-element for the Navier- 
Stokes equations. Since we neither solve the ODEs nor use numerical quadrature, our 
scheme can be precisely implemented to realize the theoretical results. It is, therefore, 
a genuinely stable Lagrange-Galerkin scheme. Our convergence results are best pos¬ 
sible for the velocity and pressure in £ 2 (L 2 )-norm for both elements as well 

as for the velocity in f°°(L 2 )-norm in the Pi+/Pi-finite element. 

The contents of this paper are as follows. In the next section we describe the 
Navier-Stokes problem and some preparation. In Section [3] after recalling the con¬ 
ventional Lagrange-Galerkin scheme, we present our Lagrange-Galerkin scheme 
with a locally linearized velocity. In Section [4] we show convergence results, which 
are proved in Section [5] In Section [6] we show some numerical results, which re¬ 
flect the theoretical convergence orders and the robustness of the scheme for high 
Reynolds number problems. In Section [7] we give the conclusions. 
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2 Preliminaries 


We state the problem and prepare the notation used throughout this paper. 

Let £2 be a polygonal or polyhedral domain of Rl 1 id = 2.3) and T > 0 a time. We 
use the Sobolev spaces LP{£2) with the norm ||-||q , W s ' p (£2) and Wq P {£2) with the 
norm || • || and the semi-norm | - | for p and a positive integer s. When p = 

2, we write If'(£2 ) = W s,1 {£2) simply and drop the subscript 2 in the corresponding 
norms. For the vector-valued function w £ W 1 “f £2 )' / we define the semi-norm \w\ i„ 
by 



o,~ 


Lq(£2 ) is the subspace of L?{£2) with the zero mean. The parenthesis (•, •) shows the 


L 2 (£2 ) ! -inner product for i = 1, d or d / d. For w £ L~(£2 ) a 
norm 


stands for the dual 


sup 

ve^(i2) rf \{0} 


(w,v) 


V 


Fora Sobolev space X(£2) we use the abbreviations H m (X) ll" , ((),T\X(£2)) and 

C(X)=C([0,r];X(fl)). 


\hdi) = {/ e HJ(t h t 2 -,H m - j (£2) d yj = 0,... ,m, \\f\\ Z m {tut2 ) < »} . 

1/2 

\\z m (t u t 2 )~ 




7=0 


and denote Z m (0, T) by Z m . 

We consider the Navier-Stokes equations: find [u,p ): £2 x (0, T) —> x K. such 

that 


— — vAu+Vp = f 

V • u = 0 

u = 0 


in £2 x (0 ,T), 

in £2 x (0, T), 
on d£2 x (0, T), 
in £2 at t = 0, 


( 1 ) 


where d£2 is the boundary of £2, + (« ■ V)n is the material derivative and 

V > 0 is a viscosity. Functions / £ C(L 2 ) and m° : £2 —► R. d are given. 

We define the bilinear forms a on Hq ( £2) d x H { j (£2) d and b on Hq (£2) d x Lq(£2) 
by 

a(u,v) = v(Vu,Vv), b(v,q) = —(V • v,q). 

Then, we can write the weak form of 0 as follows: find ( u,p): (0, T) —> Hq (£2) d x 
Lq(£2) such that for t £ (0, T), 


Du 

~Dt 



+ a(u(t),v)+b{v,p(t)) = (/(0,v), 


\/veH^{£2) d , (2a) 
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b(u(t),q) = 0, Wq€Ll(£2), (2b) 

with u( 0) = u°. 

Let u be smooth. The characteristic curve X(t',x,s ) is defined by the solution of 
the system of the ordinary differential equations, 

— (t\x,s) = u(X(t\x,s),t), t<s, (3a) 

X(s;x,s)=x. (3b) 

Then, we can write the material derivative term (4 +u ■ V)h as follows: 

+ («• v )«) (*0)v) = 

Let At > 0 be a time increment. For w : —v R. d we define the mapping X\ (vv) : 

Q -v R d by 

(Aj(w))(x) =x — w(x)At. (4) 


Remark 1 The image ofxby X\ (u(-,t)) is nothing but the approximate value of X(t — 
At\x,t) obtained by solving Q by the backward Euler method. 

Let Nj = [T /4fJ, t n = nZif and y/" = y f(-,t n ) for a function y/ defined in Q x 
(0,T). For a set of functions yr = {yr”}^ 0 and a Sobolev space X({2), two norms 

ll'llr(x) and INI flfatoX) are defined by 

\\Y\\r(x) = max {llV'lx(fl) ;n = 0,...,N T } , 

/ "2 „ \ ! / 2 

WvWp (m,n 2 X)^ 12 i 

\ n=n\ J 


and ||yr||^ (1 w r x) i s denoted by || y/j| f 2 / X j. The backward difference operator D At is 
defined by 


D A tV' = 


y/" - y f 


n -1 


At 


Let ,5* be a triangulation of Q and h = max^ ( -^ di am (AT) the maximum element 
size. Throughout this paper we consider a regular family of triangulations {<5*}/qo- 
Let V - /, x Qi, C Hq (£2) d xLq(Q) be the P 2 /P 1 - or Pi +/Pi-finite element space, which 
is called Hood-Taylor element or MINI element ||9l[2l. Let 


IlP :C{Q) d nH&(£2) d ->V h 


be the Lagrange interpolation operator to the Pj-finite element space. Let (vv*, r*) = 
r) G V* x Qi, be the Stokes projection of ( w,r) E Hq (£2) d xLq(Q) defined by 


a(wh,Vh)+b(v h ,r h ) = a(w,v h ) +b(v h ,r), Vv* £ V*, (5a) 

b(w h ,qh ) = b(w,q h ), Vg* € Qh- (5b) 


We denote by r)) 1 the first component vv* of IT* (w, r). 

The symbol o stands for the composition of functions, e.g., (go f)(x) = g(f(x)). 
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3 A Lagrange-Galerkin scheme with a locally linearized velocity 


The conventional Lagrange-Galerkin scheme, which we call Scheme LG, is de¬ 
scribed as follows. 

Scheme LG Let u Q h = (JT^(m°,0))i. Find C Vh x Qh such that 


I, 11 1, °Xi { u i, ) 
At 


Vhj + a{u n h ,v h ) +b(v h ,p n h ) = (/>/,), Vv* e Vh, 
b[u h ■ cjh ) 0, X tjfi € Qh, 


for n = 1,... ,Nt- 


Remark 2 Siili ED used the exact solution Xj‘ 1 of the system of ordinary differen¬ 
tial equations, 


dx 

-A—(t-xr)=u n h -\x;;-\t-x,n,t), t n - x <t<t n , (7a) 

X"-\t n ;x,t n )=x (7b) 

instead of X\ (w^ _ 1 ). 

By a similar way to m combined with Q, error estimates 

IK-«llr(ffi)> \\Ph~p\\p(i?) <c(h k + At), (8a) 

\\u h -u\\ r ([ 2 ) <c(h k+1 +At), (8b) 


can be proved, where k = 2 for P 2 /Pi-element and k = 1 for Pi+/Pi-element. In the 
estimate above, the composite function term (h? -1 0 X 1 v*) is assumed to be 

exactly integrated. 

Although the function n^ _1 is a polynomial on each element K, the compos¬ 
ite function u n h 1 o X\ 1 ) is not a polynomial on X in general since the image 
X\ («/ 1 ) of an element K may spread over plural elements. Hence, it is hard to cal¬ 
culate the composite function term (w? -1 0X1 (m? _ 1 ). V/,) exactly. In practice, the fol¬ 
lowing numerical quadrature has been used. Let g : K —> ffi. be a continuous function. 
A numerical quadrature //,[#; K] of f K g dx is dehned by 


Ih[g',K\ = meas(X) £ w,-g(a,-), 
i= 1 


where is the number of quadrature points and (w,,a,j £ lx /i is a pair of the 
weight and the point for i = 1,... ,N q . We call the practical scheme using numerical 
quadrature Scheme LG'. 
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Scheme LG' Let njj = (LZ|(h°,0))i. Find {( u 'ivPh)}n=i c ^ x Gfi suc h that 

KeSh 

+a(u n h ,v h ) +b(v hl p" l ) = (f n ,v h ), Vv h G V/„ 
b(u n h ,q h ) = 0 , VqhtQh, 


for n = 1,..., Nj. 

For convection-diffusion equations it has been reported that numerical quadrature 
causes the instability . For the Navier-Stokes equations we present 

numerical results showing the instability of Scheme LG' in Section[6] 

We now present our Lagrange-Galerkin scheme with a locally linearized velocity. 

It is free from quadrature and exactly computable. We call it Scheme LG-LLV. 

Scheme LG-LLV Let = (JT*(k°,0))i. Find C V/ ; x Qi, such that 

/'u1-u n r l oXiln^uT 1 ) \ 

( —---^ 1 — ~,Vh) +a(u n h ,v h )+b(v h ,p n h ) = ( f'\v h ), Vv h G V h , (9a) 

b{u n h ,q h ) = 0, \/q h G Q h , (9b) 

for n = 1....Nj. 

In the above scheme the locally linearized velocity is used in place of 

the original velocity u "~ 1 . The error caused by the introduction of the approximate 
velocity is evaluated properly in Theorems[l]and[2]in the next section. The 

following proposition assures that the integration oXi(JT^m^ _1 ),v*) can be 
calculated exactly. 

Proposition 1 Let m/„ G V/, and w G W 0 1 ' oo (f2) fl '. Suppose Ofco^Mi,” < 1, where 

CC 20 is the constant defined in ( | 12a[ > below. Then, Jq{ uh°X\ (77^ 1 w)) • Vh dx is exactly 
computable. 

Outline of the proof. When u/, and v/, are scalar functions, the result on the 
exact computability has been proved in ll23l and l22l Proposition 1], Here, we do 
not repeat the proof but show only the outline. It is necessary that the inclusion 
(Xi(Il^w))(L2) C LI holds to execute the integration of ui t oX\ (Tl^w)) ■ v* over 
LI. The condition (X 2 {\At\w\ ] < 1 is sufficient for it by virtue of Lemma[7]-(i) and 

( ] 1 2a[ > below. The mapping X\ (TJ^xv) is linear on each element. When a mapping F 
is linear, we have the following general result for any two elements Ko and K\ and 
any polynomial (/)/, of any order k defined on K\. Proposition[T]is proved by applying 
the following lemma, whose proof is easy, cf. l22l Lemma 1], 

Lemma 1 Let Kq,K\ G -9;, and F : Kq —> W 1 be linear and one-to-one. Let E\ = 
KqOF- x {Ki) and meas(£’i) > 0. Then, the following hold. 

(i) E\ is a polygon (d = 2) or a polyhedron (d = 3). 
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(ii) (j> h oF\ Ei eP k (Ei), € P*(*t). 

Remark 3 In the case of d = 2, Priestley IfTTl approximated X£~ 1 (t"~ l ',x,t n ) in (J 7 J) 

by 

Xh(x) = B\X\{x) + B 2 ^ 2 (x) + B^Xt,{x), x £ Kq 

on each K$ £ where Bj {Aj}j* =1 are vertices of Kq and {A,}^ =1 

are the barycentric coordinates of Kq with respect to {A,}^ =1 . Since X/fx) is linear in 
K(), the decomposition 

/ («r® °*h) ■ v h dx = e / («r 1 °) ■ v h dx > 

Jk ° leA(K 0 ) jE ‘ 

A (Kq) = [l-K 0 n X h \Ki) ± 0} , E,= Ko nX h \K,) 

makes the exact integration possible. However, B\ = X^~ 1 (t n ~ 1 ;A,-,t n ) are the solu¬ 
tions of a system of ordinary differential equations and it is not easy to solved it 
exactly in general since u n h 1 is piecewise polynomial. In practice, some numerical 
method, e.g., Runge-Kutta method, is required, which introduces another error. 


4 Main results 

We present the main results of error estimates for Scheme LG-LLV, which are proved 
in the next section. We first state the result when the P 2 /Pi-element is employed. 

Hypothesis 1 The solution of ([7} satisfies 

u£Z 2 nH\H 3 ), p£H l {H 2 ). 

Remark 4 Hypothesis [TJimplies (u,p) £ C(H 3 x H 2 ), which yields V • u° = 0. 

Hypothesis 2 The sequence {=^}/q 0 satisfies the inverse assumption. In addition, 
for each h, iK £ has at least one vertex in £2. 

Theorem 1 Let 14 x Q/< be the P2/Pi -finite element space. Suppose Hypotheses [ 7 ] 
Then, there exist positive constants co and ho such that if h £ (0,/io] and 
At < coh d / 4 , the solution ( Uh,Ph ) = {( u h>Ph) }'n-o °f Scheme LG-LLV exists and the 
estimates 

II Uh ~ u \\p°{ H t ), \\Ph-p\\p(j2) < ci {h 2 +At) 

hold, where c\ is a positive constant independent ofh and At. 

Next, we state the result when the Pi +/Pi-element is employed. 

Hypothesis V The solution of (JTJ) satisfies 

u £Z 2 CiH i (H 2 ), p£H l {H l ). 

Remark 5 Hypothesis V implies (u,p) £ C(H 2 x H l ), which yields V • u° = 0. 
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Hypothesis 3 The Stokes problem is regular, that is, for all g £ l?(Q Y 1 the solution 
(w,r) G H { ){£2) d x Lq(Q) of the Stokes problem, 

—vAw + 'Wr = g, 

V ■ w = 0, x € £2 

belongs to H 2 (Q) d x H l (£2) and the estimate 

holds, where c is a positive constant independent of g, w and r. 

Remark 6 Hypothesis [3]holds, for example, if d = 2 and Ti is convex |9j • 

Theorem 2 Let V/, x Qi, be the Pi +/Pi -finite element space. Suppose Hypotheses 
I 1 and [ 2 ] Then, there exist positive constants co and ho such that if h € (0,/?o] and 
At < c 0 h d / 4 , the solution ( Uh,pi ,) = {( u h'Ph)}n -0 °f Scheme LG-LLVexists, and the 
estimates 

\\ u h ~ u \\f(H l ) 1 \\Ph ~ P\\e 2 (L 2 ) <C 2 (h + At) (10) 

hold, where c 2 is a positive constant independent of h and At. Moreover, under Hy¬ 
pothesis [i] the estimate 

IK-“llr(z. 2 ) <c 3 (h 2 +At) (11) 

holds, where C 3 is a positive constant independent ofh and At. 


Remark 7 The convergence proof is easily extended for any pairs satisfying the inf- 
sup condition. However, the convergence order with respect to the space discretiza¬ 
tion is bounded by Oih 1 ) caused by the locally linearized approximation of the ve¬ 
locity. In fact, in the case of the P 2 /Pi -element the estimate ( 8 b 1 with k = 2 does not 
hold in Scheme LG-LLV, cf., Example[I]in Section[ 6 ] 


5 Proofs of the main theorems 


We prove Theorem[l]in Subsections 5.2 5.4 and Theorem[2]in Subsection 5.5 


5.1 Some lemmas 

We recall some results used in proving the main theorems. For proofs of Lemmas [2}|6] 
we refer to the cited bibliography. 


Lemma 2 (Poincare’s inequality 0) There exists a positive constant c/.\ (Q. ) such 
that 

IMIo< «i |v|!, Wv£H^(n) d . 
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Lemma 3 (the Lagrange interpolation @) Suppose { -%}h\, o is a regular family of 
triangulations of £2. Let X), be the P 2 - or Pi -{--finite element space and Tlj^ be the 
Lagrange interpolation operator to the Pi -finite element space. Then, it holds that 

imS 1 ) v||o,-<||v|| 0iOOJ Vvecooy*, 

and there exist positive constants Oho > 1 , 0!21 and OC 22 such that 


Infkkoo < a 2 ok i,oo> 

Vv£ W l '°°(£l) d , 

( 12 a) 

linZ/'v-v’ll,. < a 2 i/r" s i' 2 , 

s = 0,1, Vv £ H 2 (Q) d , 

( 12 b) 

\H l) V h \\o < 0=22 ||'7, | 0 , 

Vv/, £ Xh. 

( 12 c) 


Remark 8 The inequality ( | 12c[ > holds since X/, is finite-dimensional. If we replace 
17,' 11 by the Clement interpolation operator [[7j, this inequality holds for all v £ 
L 2 (Q) d . 

Lemma 4 (the inverse inequality fl6 ffl9) ) Suppose {■%}/,10 satisfies the inverse as¬ 
sumption. Let X/ , be the P 2 - or P\+-finite element space. Then, there exist positive 
constants 0:30 and 0:31 such that 

11V/, 11 0 ,o= < 0=30 h ~ d/6 11 V h 11 1 , Vv,, £ X h , 

|v/,|i j00 < a-i\ir dl2 |vv,|j, Vv,, £ X h . 

Lemma 5 (the inf-sup condition 02.3 .24) ) Suppose Hypothesis [2] Let V/, x Q/, C 

H^n) d xL 2 (f2) be the P 2 /P 1 - or Pi+ /P\-finite element space. Then, there exists a 
positive constant CC 4 independent ofh such that 


inf sup 

^*e< 2 /A{ 0 } v ft ev*\{ 0 } 


b{vh,qh) 

Ik/,11 1 hh I), 


> 0 : 4 . 


Lemma 6 (©) (i) Suppose Hypothesis^and that V/ ; x Qi, C Hq (Q) d x L^(Q) is the 
P 2 /P 1 - or P 1 +/P {-finite element space. Let (whHh) be the Stokes projection of (w,r) 
defined in ( 0 - Then, there exists a positive constant 0:50 independent ofh such that 

IIH, - HI !, Ik/! - Ho < a 5 oh k \\ (w, r) \\ H k+i xH k, 


where k = 2 for the P 2 /P 1 -element and k = l for the P 1 + /P \-element. 

(ii) Moreover, suppose Hypothesis [i] Then, there exists a positive constant 0:51 such 
that 

IIH-HIo < U5\h k+x \\(w,r)\\ Hk+ i xH X 

where k = 2 for the P 2 /P 1 -element and k = l for the P 1 +/P 1 -element. 

Lemma 7 (i) Let w £ W 0 1 '°°(^2) £i, and X\(w) be the mapping defined in Then, 
under the condition zif|Hi „ < 1, X\(w) : Q. —> Q. is bijective. 

(ii) Furthermore, under the condition At |w|j ^ < 1/4, the estimate 


- < det 
2 “ 


^ dXi(nQ ^ 



holds, where &Qt(dX\{w)/dx) is the Jacobian. 
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Proof The former is proved in lTT 8 l Proposition 1], We prove the latter only in the 
case d = 3 since the proof in d = 2 is much easier. Let I be the 3x3 identity matrix, 
A = (ajj) and aj = (aij,a 2 ;,£* 3 ;) r , where a,j = At dwj/dxj for i,j = 1,2,3. The 
notation | • | stands for the absolute value, or the Euclidean norm in R 3 or M 3x3 . From 
the condition 


\a\ - Oil T |i*2 + l fl 3| 2 )^ 2 + 1/4, 


we obtain 



Then, we have 



= 1 — (ail +£122 + «33) 

+ £111 £122 + £*22033 + 033a 11 —£*121*21 — <323£*32 — £*3l£*l3 — det(A) I 
<|ail + £*22 + £*33 | 

+ |<* 11 £*22 +£*221*33 +£*33 a ll —£* 121*21 ~ £*23l*32 ~ 1*311*131 + | det(A) | 

<V?>\a\ + |a| 2 + |det(A)| < 1/2, 

which implies the result. □ 

Lemma 8 Let 1 < q < °° 1 < p < l/p+l/p' = \ and w,- € Wg’°°(i2) d , i= 1,2. 
Under the condition At |w,-|i i0O < 1/4, it holds that, for y/ (42) d , 


\\\l/oXi(wi)~ V°Xi(w2)\\ 0 , q <2 1/(qp,) At\\wi-W2\\ 0iPq \\V\l/\\ 0 ’ qp ,, 
where X\ (•) is defined in (Q. 

Lemma[ 8 ]is a direct consequence of |l] Lemma 4.5] and Lemma[7]-(ii). 

Lemma 9 Let w € Under the condition zW|w|i j0 o < 1/4, there exists a 

positive constant a 6 such that, for y/ £ L 2 (Q) d , 


\\Y-Y°M^)\\-i <a 6 4f||w|| lot> Iloilo, 


where X\ (•) is defined in 

Lemma[9]is obtained from J 8 ] Lemma 1] and Lemma[7]-(ii). 


5.2 Estimates of under some assumptions 


Let 



( 13 ) 
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where ( u,p ) is the solution of <[TJ, («/,(/), Ph(t)) is the Stokes projection of (u(t),p(t)) 
defined in (0) and is the solution of Scheme LG-LLV at the step n. From (01, 

(0 and 0 we have the error equations in (efe'f): 


(D At e n h ,v h ) +a(e n h ,v h )+b(v h ,£%) = £(ft ! ,v ;i ), 

;=i 


K e l^h) = o, 


for n = \.... .Nr, where 


n _Du n u n -u n ~ l oX l (u n - 1 ) 


R' = 


Dt 


At 


n _u n - l ox l (nl l> u n - 1 )-u"- l ox l {u n - 1 ) 


R n i = 


n _n*-r,*- 1 ox l (n l 


R'i = 


At 

i^r 1 


Vv/, e v h , 
y<lh G Qh i 




ft = - 




iV 1 


(14a) 

(14b) 




(15) 


Lemma 10 Suppose Hypotheses^an ^0 Under the condition 

At\u n ~\ xi 4f|fT / [ 1) M ;;- 1 | li0 o<l/4, (16) 

it holds that 

lift 11o — ft II M llz 2 (r"— 1 ,r") > (17a) 

Iloilo </32||eri 0 + /3 3 /t 2 (||( M ,pr- 1 || ff 3 xH 2 + |« n - 1 | 2 ) J (17b) 

h 2 , 

Iloilo [\\( u , P)\\ H l(t n ~ 1 ,t n ;H 3 x H2 ) 

+ || u 'h 1 Ho,oJI( M ’7’)llL 2 (r n - 1 ,f";tr 3 x/f 2 ))> (17c) 

ii^iio<ftii«r 1 ii 0 ,.kr 1 ii ( i7d ) 

forn = 1,. .. ,N t , where pi = pi(\\u\\ c{w i.^), ft = ft(M C (wi,“) , CC 22 ), @3 = ft (McCw 1 .'”) > “ 21 , 0 : 22 , 0 : 50 ), 
ft = /3 4 (CC 50 ), ft = t/2 ant/ f/re notation pi (A) means that a positive constant ft de¬ 
pends on a set of parameters A. 


Proof We prove < 17a 1 . We decompose ft as follows: 


. , du n . . , __i . . . u" — u" x oX\(u n 1 ) . . 

7?t W =1 W + ft M ■ V)m (x)- — - - (x) 


Setting 
we have 


+ ( u n (x) - «" - ■ 1 (x) ) • Vm” (x) = ft'! (x) + ft 2 (x). 


y(x,s) = x+ (s — \)At u n 1 (x), t(s)=t n l +sAt, 


/ l -u n - 1 oXi(u n ~ l ) 1 


4r 


= At [“W+)’ f ( s ))]l 


s=0’ 
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which implies that 


du ' 1 


dt 

i' 1 

= At ds 


du 


R"n = ^r + («" • v)«"- /, + 





= At Jo Sl j(“" 1 (‘)’ V + ^) “| CK-.-siy^OKsi- 


Hence, we have 


*ii Ho < At / si 


«" (')' v + 57 ) u f M-^OtOi)) 


t/si 




where we have used the transformation of independent variables from x to y and si to 
t and the estimate | det(<9x/cFy)| < 2 by virtue of Lemma[7]-(ii). It is easy to show 

du 

L 2 (t"-\t"-L 2 ) 


||*? 2 || 0 <VS|h"|i,. 


From the trian gle in equality we get ( |_F7a| >. 

We prove (17b i. Using Lemma |8| with q = 2, p = 1, p’ = °°, wi = 


W 2 = u n 1 and t \f = u n 1 , we have 


Iloilo <\u n ~ l 




From Lemmas|3]and[6]-(i) we evaluate the first term as follows: 

linfuT 1 -niV-'Ho = UnW^r 1 -n^V-^Ho 
<oe 22 ||mS _1 

<«22(||«r 1 -^ 1 || 0 + ll ^ 1 IIo + IIm ”- 1 -n'V'io 

<a22(\\el- l \\o + a 5 oh 2 \\(u 1 p) n - l \\ H i xH 2 + a 2 ih 2 \u n - l \ 2 ). 

The second term is evaluated as follows: 

\\n^u n - l -u n -%<a 2 xh 2 \u n - l \ 2 . 

Thus, we have 

ll*2l|0 < |‘ | ! >00 {0^22(11*|o + a 50h 2 \\(u,P)"~ 1 \\h 3 xH 2 ) 


+ a 2l {l + a 22 )h 2 \u n 1 1 2 ^, 


( 18 ) 
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which implies ( 17b| >. 
We prove (| 1 7c[>. Let 


y(x) = x+(s — l)AtTI^u" l {x), t{s) = t n l +sAt. 


Since R'\ is rewritten as 


R1 = / o 1 {^ + (n( 1 ) M "- 1 (.)-V)r,J (y(., s),t(s))ds 

have, by using the change of the variable and Lemma[7]-(ii), 

1 dt] 


o< 


< 


dt 


(y(-,s),t(s))ds 




( 

<9 77 


dt 


\J~At V 
V2a 5 oh 2 
\f~Kt 

which implies 


< 


+ ||Lr( 1) «”- 1 llo,=c||V77ll i 2 (; „-, ir „ ;i 2 ) 

L 2 (t n ~ l 

(\\( u iP)\\H'(t"- l ,t n :H 3 xH 2 ) + hi, * ||o,oo II ( u i P)\\L 2 (t n ~ l ,t":H 3 xH 2 )^ 


ich implies ( |17c[). 

The inequality (|17d i is obtained from Lemma[ 8 ]with q = 2, 


W 2 = n ( h 1 ] u n h 1 and l// = e n h 1 


= 2, p = °°, p' = 1, wi = 0, 
□ 


Lemma 11 Suppose Hypotheses [7]aH£/[2] Let n € {1, • • • . Nj } be any integer and let 
ul l € Vh be known. Suppose that 1 satisfies 

b(u n h ~\q h ) = 0, V®e6, (19) 

Under the condition © there exists a solution of (|9| and it holds that 

\\D At e'lf 0 +D A ,(v\el\\) 

</fc 1 (|K- 1 || 0 ,»)v|^- 1 |? + fe(||«2- 1 || 0 , oo ){^ll«ll| a( ^i^ ) 

+ ^( U,P ^H 1 (t"~\t"-Jl 3 xH 2 ) +h4 (\\( U ’P)" 'll H 2 xH 2 + K * 1 2 ) }’ 


where e n h is defined in (13), and /3t 1 (<5) and are the functions defined in (21) 

below. 

Proof Since it holds that At\n^' > u n h ~ l \\^ <1/4, the mapping X\ (JT^kJJ -1 ) : £2 — ► 
£2 is bijective from Lemma 7j-(i). Hence, there exists a solution (u'^pf) of ([9|. Sub¬ 
stituting v/, = D At e n h in (14a >, we have 


\\D A ,el\\l+D At ( V - HVeJIlg) +b(D At ele n h ) < £(/?",L>^)- (20) 

y 1=1 
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From ( fl9| ) and ([9]» the term b(D At e'l l .e" l ) of the left-hand side vanishes. Using Schwarz’ 
and Young’s inequalities and Lemma[lO] we have 


ll^r^Ho +D At < 2|/3| At 

+(/3 2 || e r 1 ii 0 +/33/t 2 (ii(«,pr l n ff 3 Xff2 +| M "- i | 2 )) 2 


+Ar— 
4 At 


P)\\H l (t n ~ l ,t n ;H 3 xH 2 ) + H 1 || 0 ,oo ll(“i/ 7 )llz, 2 (f"- 1 ,r";/t 3 x// 2 )) 

+Ps\K%,K'\1} + 1\\Da,4\\1, 


which implies that 

£ , 4f(v| e /!ll)+ ll^ztr^llo 

<j8h (ikr'llo+ll^r 1 iio,~ kr 1 1?) 

/i 4 2 1 2 2 

+ ~^(\\(. u iP)\\H l (t"- l ,t n -,H 3 xH 2 ) + \\ U h llo.oo ll( M )/ , )llL 2 (f'>- 1 ,/";H 3 X /f 2 )) 

+h 4 (\\( u , P y-% 3xH2 + \ u ''-'Q}, 


where pi \ and p \2 are constants depending only on pi.... .p^. Using Poincare’s in¬ 
equality ||e^ -1 ||o < OCi|e^ _I |i and defining the functions j3 2 i and j3 22 by 

/ 3 2I (^) = ^ 1 (a 1 2 + ^) ! / 3 22 (^) = / 3 12 (H- 4 2 ), (21) 

we have the conclusion. □ 


5.3 Definitions of constants c*, co and Iiq 
We first define constants J3|j and j3* 2 by 

p 2 \ = /32l(||w||c(L“) + 1), P 22 = ^22(||m||c(L“) + !)■ 

We define two positive constants c* and co by 

c* ={v~*(l + af)exp(/3|ir)/3 22 } 1/2 max{||M|| z 2 , (||(M,p)||^i (ff 3 xH 2 ) 

+ ^ ^11 ( M ’P) llc(/T 3 xH 2 ) + Hc(ff 2 )) + Va 5o|p 0 |2) } 

and 

C0 =I / 1 

C ° 4 V OhoO^lC* ’ 

Let a positive constant /zo be small enough to satisfy that 

«3ofy) J/,fi 4“ Ct5o/zo|| iPiP) Wc(H 3 xH 2 ) + a 21 Mc(tf 2 )) 
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-\-aj,QC it CQh cl J n < 1 , 


C()|«3i/!q C> ^ (ct,h Q + a 5 oh 0 \\(u i p)\\ C ( H 3 xH 2' j +«21 |«lc(i/2 

+Ci2ohJ — 

which are possible since all the powers of /zq are positive. 


16CC20 


(23a) 


(23b) 


5.4 Induction 

For n = 0,.... A'y we define the property Pin) by 

(a) v \e n h \ 1 + || D A ,e h \\ 2 e2{l n . Ll) < exp \\u\\ 2 Z 2 {t o r) 

+ h4 (\\( U ’P)\\ 2 H'(P,t"-,H 3 xH 2 ) + ll( M !F)ll^(0 in _l; ff 3 X H2) + \ U \fi(0,n-UH 2 )) 

+ V \ e h\l}' 

(b) ||<|| 0jOO < I|m|Ic(L“) + 1 - 

(c) At\n^ ] u n h \xr < 1/4. 

Proof of Theorem^ We first prove that P (n) holds for n = 0,... ,Nt by induction. 
When n = 0, the property P(0)-(a) obviously holds with the equality. The properties 
P(0)-(b) and (c) are proved in similar ways to and easier than P(«)-(b) and (c) below, 
we omit the proofs. 

Let n £ {1, • • • . Nr} be any integer. Supposing that P (k), k = \.... ,n — 1, holds 
true, we prove that P (n) holds. We now apply Lemma [TT| The condition ( fl9| is sat¬ 
isfied trivially when n> 2. When n = 1, from the choice of |5jl and Remark]^] we 
have 

b{eh,qh)=b(u° h ,q h )-b{u° h ,q h )= 0-0 = 0, \/q h £Q h ■ (24) 

We consider the condition ( fl6| . The former condition follows from At < C()h d ' 4 and 
( |23b[ ) by the inequality 

At\n n 1 1 1,00 < cqIiJ \u\c(w l ’°°) — i6 a20 — 

and the latter condition AtlPI ^< 1/4 follows from P(« — l)-(c). Hence, 
there exists a solution (u'} v p",) at the step n. 

We begin the proof of P(n)-(a). By putting 

2 11_ 11 2 

x n=y\ e h\\i yn = H-CLne/,|| 0 , 

bn =At |Mlz 2 (r"- 1 1 f' 1 ) 

+ b 4 (iill(“’P)ll//l(f"-l,;";i/3x// 2 ) + ll( M ’P)” 1 \\h 3 xH 2 + l“" * I 2 ) ’ 

P(n)-(a) is rewritten as 

„ +At < exp(/3 2 *| nAt)fi 22 (x 0 + At £ b\ 
i= 1 V ;=i / 


x, 


(25) 
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On the other hand. Lemma [TT| implies that 

x n + Aty n < (1 + ^21 At)x n -i + P 22 Atb n , 

where we have used the inequalities /32i(||t<^ _1 |L oo ) < i = 1,2, obtained from 
P(n — l)-(b). Using the inequalities 1 < 1 +x < exp(x) for x > 0 and P {n — l)-(a) 
rewritten by ( |25) l, we have 

n n— 1 

x n +AtY j y i =Xn+Aty„+AtY j yi 

i= 1 i=l 

n— 1 

<(1 + $ 2 \At)x n ~i + ^ 22 Atb n + 

i=l 

<(1 + /3 2 |4f)exp(/3 2 i (n — 1 )4?)/3 22 ^ bi'j + /3 22 Atb n 

< exp(/3 2 Vi4?)/322E^ > 

which is nothing but P(u)-(a). 

Since is the first component of J7 *(m°, 0), we have 

el = u° h -^ h = (77^(0 = ((n* -/)(0, -p°))t, 

which implies < aso/i 2 |p°| 2 . From P(0)-(a) and the definition of c*, we have 

\\e n h \\ X <c*{h 2 + At). (26) 

P(n)-(b) is proved as follows: 


Kilo,- 

<M- n h )un \\o,~ + \\ n h ) u n \\o,oo 

<a 30 /t -d / 6 |K -JT^VlIi + ||m"|| 0 oo (by Lemmas [4] and [3]) 

<0E 30 /*- rf/6 (||«S - ^11! -Hl^ - «"ll t + ll«" - -rri 15 **" Hi) ll«"llo,o- 

<a 30 h~ d/6 (c,(/t 2 + 40 + a5o/t 2 ||(M,p)l/f3 xH 2 + a2i/tK| 2 ) + ll«lo,oo 

(by ( p6| ), Lemma|6]-(i) and Lemma[3]) 

<a 30 h l ~ d/6 (cJi + a 50 h\\(u,p) n \\ H3> , H 2 + a2i |m"| 2 ) + a3oc*c 0 // /12 + ||m"|| 0oo 

(since At < coh d / 4 ) 

<1 + ||m||c(l“) ■ (since h < ho and by ( |23a| >) 

We prove P(n)-(c). We can estimate |«J'| 1 

Kli,«A 

<{\u n h -n^u n \^ + \n^u n \^)At 


At as follows: 
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+ OC20 \u n \ la> }At (by Lemmas [ 4 ] and [ 3 J 

<{a 31 h~P 2 (\u n h - SJIj + \tT h - u n |j + |m" - JT A (1) « n |i) + «20 
<|«31^ _d/2 (c*(/! 2 +4f) + «50/2 2 ||(M,p)''||//3 xi/ 2+a2l/!|i("l2) 

+ CC20 \u’\ oo|^ f (by ( p6| ), Lemma[6]-(i) and Lemma[ 3 ]) 

<co|a3i/! 1_d/4 (cji + a 5 Q h\\(u, p) n \\ H i xH 2 +a 2 \ \u n \ 2 ) + a 2 oh d/4 \u n \ llx j 

(since At < C(,h d / 4 ) 


O£31 c * c 0 


<- 


1 


1 


1 60:20 160:20 4 a 2 o 


(since h < ho, and by ( 23 b[> and the definition of cq) 


From this estimate and the definition of Ohio we have At\Ilj - 1/4- 

Thus, we have proved that P(n) holds for n = 0, • • • ,Nt- 
From P(n)-(a), n = 0,... ,Nt , we obtain 

Ik/ill^^ 1 ) > ||-D4f^/i||^(i jAfj , ;i 2) <c*(/z +At). 

Using the triangle inequality \\u h -u\\ r{Hl) < \\e h \\ r{H \ ) + ||t? we get 

< c\(h 2 +At). 

We now prove the estimate on the pressure. We can evaluate e", as follows: 

(by Lemma [5]) 


|e;il„<4 sup 


°kt v/,ev* ||v*||i 


1 

= — sup 


1 


vfjGV/, Ir/ill 1 \i=i 

4 


E ’ v * )-ipA t e n h ,v h )-a ( e n h , v h ) 


(by ( 14 a 1) 


<^(eii^IIo+II^|| 0 +vk;i 1 ) 

< c ^||r>4f^|| 0 + V |eJJ|j + ||e^ 1 111 + \\u\\ z 2( t n-ljn) +/t 2 ||(«,p)" 1 ||^3 x ^2 

h 2 9 \ 

+ ^^\\( u ’P)\\ Hl (t n - 1 ,t n fi 3 >< H2 )+ h 2 \u n 1 1 2 J , 

(by Lemma 10 and P(n — l)-(b)) 

which implies that, from ( [26] ), 

||e/;||^(L 2 ) < c(||r ) 4;e/ J ||^2( i 2) + h~ + At) < c(/r + At), 
where c is a positive constant independent of h and At. Using the triangle inequality 
IIP/! ^Pll^(L 2 ) < ll e A|l^(L2) + l|P-Pft||^2( L 2), 
we obtain \\p h -p\\pr L 2 ) < c\(h 2 +At). □ 
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5.5 Proof of Theorem [2] 

In this subsection we prove the result on the Pi+/Pi-element. At first we replace the 
estimates of R" and R" in Lemma [T0| 

Lemma 10' Suppose Hypotheses 1' and[2] Under the condition ( | 1 6[ i it holds that 

Iloilo < ft ||er 1 || 0 + j3 3 (/t||(«,pr- 1 L2 x ^ +h 2 |m" _1 | 2 ), 

Iloilo — ft^^ \\H l (t n ~ 1 ,t";H 2 xH 1 ) + II u h 1 Ilo,oc II ( M > P) \\L 2 (t'‘- 1 ,t";H 2 xH' )) 

for n = 1 ,... ,Nt- 

The proof is similar to Lemma [T0| by replacing the order k = 2 by k = 1 in Lemma 

B®. 

I 1 roof of Theorem [2] We only show the outline of the proof for the existence of (uh-Ph) 
and the inequality (|1()|> since the proof is similar to that of Theorem[l] We replace the 
definition of c* by 

c* ={v _1 (l+ai 2 )exp(^| 1 r)^ 2 * 2 } 1/2 max| \\u\y,(\\(u,p)\\ 2 Hl{H 2 xHi) 

+ T (ji\i u ,P)\\c(H 2 xH l ) + MC(// 2 )) + Va 5ol/ ;,0 |l) }) 

redefine co by ( (22) with the new c*, and replace the condition ( |23) on ho by 

^30^0 ^ “I - t^501| ( W 5/?) llc(// 2 x// 1 ) "b ^21 Mc(// 2 )^ ~b ^30 c * c 0^(/ — (27a) 

co { « 3 1 / 2 o " //4 (c* + a 5 0 II (m , ft ||c( H 2 xffl) + «21 | m| c(H 2 } ) 

+ « 2 ofy/ l M lc(tv 1 '°°)} — ■ (27b) 

We also replace P(«)-(a) by 

V l e /!ll + ll^ 4 ' e/! ll«2(l,,i;L 2 ) — ex P(ftl ,l ^0ft:2 IMIz 2 (f°,f") 

+ ll~ (\\(u,p)\\ 2 H i^O :t n . H 2 x H l) + \\( U ’P)\\fi( 0 ,n-l;H 2 xH l ) + \ U \fi( 0 ,n-VJl 2 )) + V Kll}' 
P(n)-(a) implies the estimate 


Kilt <c*(h + At). 


(28) 


The choice ( |27| ) is sufficient to derive P(n)-(b) and (c). Hence, the existence of the 
solution and the estimate ( [TO} are obtained similarly. 

We now prove the estimate jlT I, following fl9l except the introduction of X\ (17^ 11 1 ). 

Substituting (■ Vh,qh ) = { e h’ £ h) i* 1 * 14 1 , we have 


^DMl + ^Wel-el 


i + v|«J| ? = K«), 

1=1 


(29) 
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where R t , i = 1, • • • ,4, are defined in The term (R".e" t ) is evaluated by ( |17a| >. 
From Lemma[ 6 [(ii) we have 

-u n -'\\ 0 < a 5 ih 2 \\(u,p) n ~ l \\ H 2 X Hi- 

Using this estimate in the last line in ( [18] ), we have 

(RiA) < Iloilo Iloilo < {ft I|^" 1 |Io+^ 2 (II(“.p) b_i II^x«i + |«"" 1 | 2 )}II^IIo- 

We divide the term as follows: 

wa) =± ( r - r i" 1 , e h )- r- 1 o x,(n< 

+ j t (r - 1 o Xl (n fc (1 V- 1 ) - v" 1 ox, 

= I\ + I 2 +h- 


The first term I\ is evaluated as 


h < 


\/~Kt 


<9 77 


dt 


L 2 (t n ~ l ,t";L 2 ) 


e h\\o — ~^r^~\\( U 'P)\\H l (t n - l ,t";H 2 xH') ll e /!llo • 


By Lemma|9]the second term I 2 is evaluated 


as 


h <«6O!20 ||m|Ic(W 1 ’“) II n” lolKlll 

<a 6 a 20 \\u\\ c{wln a 5 ih 2 \\(u,p) n - l \\ H 2 xHl IKI^. 


In order to evaluate / 3 we prepare the estimate 


«30 h d/6 \ rj" 1 1 j < a^oh 1 d/6 a 50 \\(u,p) n l \\ H 2 


xH l ■ 


< 1, 


where we have used Lemma[ 6 ]-(i) and (|27[i). 

Using Lemma 8 [ w ith q = 1, p = p’ = 2, w\ = TZ^'V , n „ 

\jf = Tj n ~ 1 , Lemmaffl the above estimate and (18 1 , we can evaluate I] as follows: 


W2 = n { h l) u‘ 


and 


h<l~ t ||t7 , '- 1 ox 1 (n| 1) M"- 1 )-r7''- | ox 1 (n‘ 1) M r 1 )||o a ||^|| 0iTO 

<^a 3 o/t- d/ 1r 7 ”- 1 | I ||n(\«- 1 -niV,- 1 ||o||4llt 

<V2a22{\\e^~ 1 \\o + a5ih 2 \\(u,p) n ^ l \\ H 2 xH i + a 2 \h 2 |m” — 1 1 2 ) H^llt ■ 
In order to evaluate (R'^e 1 },) we prepare the estimate 

“ 30 h~ d/6 < awh~ d/ 6 c*{h + At) < CC 3 oc*h~ d / 6 (h + coh d/4 ) < 1 , 
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where we have used ( [28] ) and (f27}i). Using Lemma |9j the above inequality and a 
similar estimate to It, in (R",e'p), we can evaluate as follows: 


(Rlel)=-j- t (el- l -er 1 oX l (n^u n - l ),e n h ) 

<a 6 a2o||M|| C (tv^)lkr 1 llolkftlli 

+ V 2 a 22 (||cr 1 llo + «5i/t 2 ||(M,p) n “ 1 ||//2 xH i+a2i/t 2 |M"“ 1 | 2 )||^llf 


Combining (291 with these estimates and using Young’s inequality and Poincare’s 
inequality 11 ^*T[o < a i | e /i|p we have. 


Dm K\\l + v |ej*|j < fti 1 1| 0 + \TKt ||m|| Z 2 ^„-i^ 


1 ^ 

2 


,t") 


+ ^ t \\( u , p )\\H'(t"-l,t"-M 2 xH')+ h 2 U u iP) n 1 \\h 2 xH' +h 2 \u n ’| 2 ) \\e n h \\ J 


1 n— ] 

LI \ 

\u 

I2) 


2 || i 112 f 2 /t 

<V Wl t \l + p32 \\e n h ~ || 0 + fe|4r ||M||" 2 ( r „_l jf „) + — \\{u,p)\\ H l(,n-l r . H 2 xH 1 ) 

+ h 4 (\\{u,p) n 1 II^2 X ^1 + \u” 1 1 2 ) }> 

where ft 31 , ^32 and ft 33 are positive constants independent of h and At. Applying 
Gronwall’s inequality, we obtain {FT}. □ 


6 Numerical results 


We show numerical results in d = 2 for the P 2 /Pi-element. We compare the con¬ 
ventional Scheme LG'with the present Scheme LG-LLV. For the triangulation of the 
domain the FreeFem++ HD is used. In Scheme LG'we employ numerical quadrature 
of seven-point formula of degree five ED. The relative error Ex is defined by 


ExW 


1 QhWx 

IWWx 


for (j) = u in X = £°°(Hq ) and £ °°(L 2 ), and for (j) = p in X = £ 2 (L 2 ). 

Example 1 In ([lj, let £2 = (0, l) 2 , T = 1. We consider the two cases, v = 10~ 2 and 
10 4 . The functions / and u° are defined so that the exact solution is 


mi {x,t) = 0(xi,x 2 ,O, 
u 2 {x,t) = -0(x 2 ,xi,f), 
p(x,t) = sin(7r(xi + 2 x 2 ) + 1 +t), 

where <j)(a,b,t) = — sin(7ra) 2 sin(7r/?){sin(7r(a-|-f))-F3sin(7r(a-|-2Z? + f))}. 
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Fig. 1 The triangulation of Cl for N — 16 

Let N be the division number of each side of Q. We set h = 1 /N. Figure fl] shows 
the triangulation of Q for N = 16. The time increment At is set to be At =Jr (N = 
16,23,32,45 and 64) or At = /r 3 (N = 16,19,23,27 and 32) so that we can observe 
the convergence behavior of order h 2 or The purpose of the choice At = 0(h 2 ) 
or O(Jr') is to examine the theoretical convergence order, but it is not based on the 
stability condition, which is much weaker as shown in Theorem[l] 


Table 1 Symbols used in Figs. 


[^Jandjs] and TablesJ^Jandjs] 


0 

u 

P 

u 

u 

X 

r(HD 

f(L 2 ) 

r(L 2 ) 

r(L 2 ) 

At 

h 2 

li 2 

h 2 

h 3 

Scheme LG' 

• 

■ 

▲ 

▼ 

Scheme LG-LLV 

O 

□ 

A 

V 


Table 2 The values of relative errors and orders in Fig. |2] by Scheme LG'(top) and Scheme LG- 
LLV(bottom) 


N 

• 

order 

■ 

order 

▲ 

order 

N 

▼ 

order 

16 

8.55e-2 


1.63e-l 


7.77e-2 


16 

6.45e-3 


23 

4.34e-2 

1.87 

8.40e-2 

1.82 

4.03e-2 

1.81 

19 

3.73e-3 

3.19 

32 

2.30e-2 

1.93 

4.52e-2 

1.88 

2.17e-2 

1.87 

23 

2.10e-3 

3.02 

45 

1.20e-2 

1.90 

2.34e-2 

1.92 

1.13e-2 

1.93 

27 

1.29e-3 

3.02 

64 

6.02e-3 

1.97 

1.18e-2 

1.96 

5.64e-3 

1.96 

32 

7.57e-4 

3.15 


N 

O 

order 

□ 

order 

A 

order 

N 

V 

order 

16 

8.97e-2 


1.93e-l 


7.84e-2 


16 

1.48e-2 


23 

4.62e-2 

1.83 

1.03e-l 

1.73 

4.10e-2 

1.78 

19 

9.19e-3 

2.78 

32 

2.46e-2 

1.92 

5.44e-2 

1.92 

2.25e-2 

1.82 

23 

6.04e-3 

2.19 

45 

1.29e-2 

1.90 

2.84e-2 

1.91 

1.17e-2 

1.93 

27 

3.83e-3 

2.85 

64 

6.39e-3 

1.99 

1.41e-2 

1.97 

5.81e-3 

1.98 

32 

2.72e-3 

2.01 


Table|T]shows the symbols used in the graphs and tables. Since every graph of the 
relative error Ex versus h is depicted in the logarithmic scale, the slope corresponds 
to the convergence order. Figure |2l shows the graphs of (p) an d 
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h h 

Fig. 2 Relative errors (p), ( M ) — h 2 (left) and Eg^^i^ ( u ) with At = /z 3 

(right) in the case of v = 10~ 2 in Example]!] 


Table 3 The values of errors and orders of the graphs in Fig. p] with At — h 2 (top) and At = /z 3 (bottom) 


N 

• 

order 

■ 

order 

▲ 

order 

N 

▼ 

order 

16 

1.91e+0 


2.14e-l 


1.93e-l 


16 

1.55e-l 


23 

1.34e+0 

0.97 

8.97e-2 

2.39 

8.81e-2 

2.16 

19 

6.64e-2 

4.92 

32 

9.42e+0 

-5.90 

3.48e-l 

-4.11 

5.28e-l 

-5.43 

23 

3.65e-2 

3.14 

45 

4. lOe+1 

-4.31 

1.28e+0 

-3.81 

1.46e+0 

-2.98 

27 

1.92e-2 

4.01 

64 

B.B2e+l 

-2.18 

2.77e+0 

-2.20 

2.02e+0 

-0.93 

32 

1.02e-2 

3.71 


N 

O 

order 

□ 

order 

A 

order 

N 

V 

order 

16 

6.72e-l 


2.65e-l 


2.09e-l 


16 

2.47e-l 


23 

3.91e-l 

1.50 

1.36e-l 

1.83 

9.88e-2 

2.07 

19 

1.05e-l 

4.96 

32 

1.85e-l 

2.26 

6.98e-2 

2.02 

4.18e-2 

2.60 

23 

8.80e-2 

0.94 

45 

1.27e-l 

1.10 

3.73e-2 

1.84 

2.12e-2 

1.99 

27 

6.18e-2 

2.20 

64 

7.21e-2 

1.61 

1.83e-2 

2.03 

9.78e-3 

2.20 

32 

2.97e-2 

3.29 


ZW 2 )(ii) versus h in the case of v = 10 2 . Their values and convergence orders 
are listed in Table |5j When At = h 2 , the convergence orders of E e ^^ H i^(u) (•, O), 
Ef 2 , L 2 ) {p) (■, □) and E^, L 2 -, (u) (A, A) are almost 2 in both schemes. When At = /j 3 , 
the order of E e «,/ L 2 \(u) is almost 3 in Scheme LG' (T) and 2 in Scheme LG-LLV (V). 
They reflect the theoretical results. 

We consider a higher Reynolds number case. Figure [3] shows the graphs in the 
case of v = 10 4 and their values are listed in Tabled When At = h 2 , all errors 
increase abnormally at IV = 32,45 and 64 in Scheme LG' (•, ■. A) while the con- 
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64 32 16 

h 

Fig. 3 Relative errors E £o0 ^ (w), E £2 ( L i ) (p), 
(right) in the case of v = 10~ 4 in Example^ 



64 32 16 

h 

(L 2 ) ( M ) At — h 2 (left) and E £OQ ^ L 2 ^ (u) with At — /z 3 



Fig. 4 The domain Q. and the boundary condition (left) and the triangulation of El (right) in 


Example 


s 


vergence is observed in Scheme LG-LLV (O, □, A) but the order of (O) 

is less than 2. In order to obtain the theoretical convergence order 0(h 2 ) in Scheme 
LG-LLV, it seems that finer meshes will be necessary. When At = h \ the order of 
A”(l 2 )( m ) i s mo re than 3 in Scheme LG' (T) while it is less than 3 between /V = 19 
and 23, and N = 23 and 27 in Scheme LG-LLV (V). 

We now consider a cavity problem to see that Scheme LG-LLV is robust for high 
Reynolds number while Scheme LG' is not. This problem is not a homogeneous 














24 


Masahisa Tabata, Shinya Uchiumi 



Fig. 7 



Dirichlet boundary problem, but it is often used as a benchmark problem. In order to 
assure the existence of the solution we deal with a regularized cavity problem, where 
the prescribed velocity is continuous on the boundary. 

Example 2 Let £2 = (0, l) 2 , / = 0, w° = 0. We consider the two cases, v = 10 -4 and 
1CL 5 . The boundary condition is described in Fig.[4](left), where g\ = 4x\ (1 — xj). 

Figure [4] (right) shows the triangulation of £1. Figures [5] and [6] show the stere¬ 
ographs of the solution u n h at t n = 8 in the subdomain (0.3,0.7) x (0.8,1.0) by Scheme 
LG' and Scheme LG-LLV, respectively, when v = 10 4 . Neither solution is oscillat¬ 
ing although Uhl of Scheme LG' takes larger values than that of Scheme LG-LLV. 
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Figures [7] and [8] show the stereographs of the solution u n h at t" = 8 in the subdo¬ 
main (0.3,0.7) x (0.8,1.0) by Scheme LG' and Scheme LG-LLV, respectively, when 
V = 1 0 5 . While oscillation is observed for both components of the solution by 
Scheme LG' in Figure [7] we can see that the solution by Scheme LG-LLV is solved 
without any oscillation in Figure [8] 


7 Conclusions 

We have present a Lagrange-Galerkin scheme free from numerical quadrature for 
the Navier-Stokes equations. By virtue of the introduction of a locally linearized 
velocity, the scheme can be implemented exactly and the theoretical stability and the 
convergence results are assured for practical numerical solutions. We have shown 
optimal error estimates in £°°(H l ) x £ 2 (L 2 )-norm for the velocity and pressure in the 
case of P 2 /P 1 - and Pi+/Pi-finite elements. Numerical results have reflected these 
estimates and the robustness of the scheme for high-Reynolds number problems. 
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